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Abstract 

We present a class of optimum ground states for spin-| models on the 
Cayley tree with coordination number 3. The interaction is restricted to 
nearest neighbours and contains 5 continuous parameters. For all val- 
ues of these parameters the Hamiltonian has parity invariance, spin-flip 
invariance, and rotational symmetry in the xy-plane of spin space. The 
global ground states are constructed in terms of a 1-parametric vertex state 
model, which is a direct generalization of the well-known matrix product 
ground state approach. By using recursion relations and the transfer ma- 
trix technique we derive exact analytical expressions for local fluctuations 
and longitudinal and transversal two-point correlation functions. 



1 Introduction 

The Cayley tree belongs to the category of pseudo-lattices Unlike regular 
lattices, which are usually defined in terms of periodic structures, the Cayley 
tree is generated by the following recursive scheme: 

1. A Cayley branch of order 1 is a single lattice site. 

2. A Cayley branch of order n is defined as a lattice site with K+l bonds, to 
which K Cayley branches of order n— 1 are attached, i.e. the branch has 
1 unconnected bond. 

3. A Cayley tree of order n is given by a central lattice site with K+l bonds, 
to which K+l branches of order n—1 are attached. 

K is called the connectivity, K + 1 is the coordination number of the Cayley 
tree. Figure [l] shows a finite Cayley tree with coordination number 3. In the 
thermodynamic limit n — > oo the Cayley tree is also known as the Bethe lattice. 

An important property of the Cayley tree is that there is exactly one path 
from a lattice site i to another lattice site j. The distance \i—j\ between i and 
j is simply the number of edges on this path. We can use the distance function 
to divide the Cayley tree into two disjoint sublattices. Denote the central site 
by iq. An arbitrary lattice site i belongs to the sublattice La if |*o — i\ is even, 
otherwise i belongs to the sublattice Cb ■ It is easy to see that every site i € Ca 
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Figure 1: A finite Cayley tree with coordination number 3. 



has only nearest neighbours in Cb and vice versa, so (Ca,£b) is a bipartite 
decomposition of the Cayley tree. Note that it is valid for any connectivity K . 

Due to the hierarchical structure of the lattice, the partition function of 
a classical statistical model on the Cayley tree can usually be calculated by 
using recursion relations, provided the local interaction has finite range and 
the number of states at each lattice site is also finite. This is a substantial 
step beyond models on the chain, as it allows to construct exactly solvable 
models with arbitrary coordination number. The most important drawback of 
the Cayley tree is that for large system sizes -^th of the lattice sites are boundary 
sitesQ. As a consequence the physics is heavily influenced by boundary effects. 
There is no canonical way to impose periodic boundary conditions. 

In this work we investigate the ground state problem of a class of quantum 
spin- 1 models on the Cayley tree with connectivity K = 2. The paper is 
organized as follows. Section |2| contains the definition of the Hamiltonian and 
discusses its parameters and symmetries. The global ground state is constructed 
explicitly in section ^ in terms of a vertex state model. Vertex state models are 
graphical realizations of so-called optimum ground states, which simultaneously 
minimize all local interaction operators. Ground state properties, i.e. single- 
spin and two-point expectation values are presented in section ^. As shown 
in appendix A the calculation of ground state expectation values leads to a 
classical vertex model on the Cayley tree, which can be solved exactly. Finally 
we summarize our results in section pi 



For small system sizes the percentage of boundary sites is even larger. 
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2 Model definition 



The model is defined on the Cayley tree with coordination number 3. A spin- 
| is located at each lattice site. These spin variables are coupled by nearest 
neighbour interaction terms hij, which are all equal, they only act on different 
pairs of lattice sites. The global Hamiltonian H is the sum of all these local 
interactions, so the system is completely homogeneous. 

The local interaction is the same as in our previous works on the hexagonal 
lattice H and the two-leg ladder Hence we shall be very brief here. For the 
construction of optimum ground states it is advantageous to write the interaction 
operator in terms of projectors onto its eigenstates 



hij = A 3 (\v 3 )(v 3 \ + \v-3)(v-3\) + 

A 2 - ff (\v^)(vr\ + \vZ a 2 )(vZ a 2 \) + 
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If we use the following notation for the canonical basis states of a single spin- 
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The parameters A3, X 2 a , Aj" 2 , A^ 7 are real and positive and the superposition 
parameter a is real, a is a discrete parameter, which can only take the values 
±1. Thus the total number of continuous parameters is 5, which includes a 
trivial scale, so there are 4 non-trivial interaction parameters. 

For all values of the parameters (Q) commutes with the pair magnetiza- 
tion operator Sf + Sj and with the parity operator Py , which interchanges the 
spins at sites i and j. Therefore the local interaction (Q) has rotational symmetry 
in the xy-plane of spin space and is parity invariant. In addition, corresponding 
eigenstates with magnetization m and — m carry the same A-coefhcient, so hij is 
also invariant under a spin-flip S z — > — S z . In particular, no external magnetic 
field is applied. 

As all A-parameters are positive, (|l|) is a positive semi-definite operator, i.e. 
all its eigenvalues are non-negative. The two-spin states (||) are the excited local 
eigenstates of h^, the remaining 9 eigenstates are local ground states, i.e. the 
corresponding eigenvalue is zero. Since the Hamiltonian H is the sum of positive 
semi-definite operators, the global ground state energy Eq is non-negative, too. 
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In the next section we shall show that Eq is in fact zero and the corresponding 
global ground state will be constructed explicitly. 

At the isotropic point, a — — -\/3 and a = — 1, the A-parameters can be 
adjusted so that hij has the form 

»«-S,.S J + ^(B ( .B J )' + i(B,.S J )' + ^. (4) 

Obviously, this operator has complete 50(3) symmetry. It simply projects onto 
all states with (S^ + Sj) 2 = 3(3 4- 1). This case has already been investigated in 
0. Its ground state is known as the valence bond solid (VBS) ground state. As 
shown in it has exponentially decaying correlation functions, no Neel order, 
and there is an energy gap between the ground state and the lowest excitations. 
This is consistent with our results presented in Section ^. 

3 Construction of the global ground state 

In this section we construct the exact ground state of the present model. It 
is an optimum ground state Q-Q, i.e. it is not only the ground state of the 
global Hamiltonian H, but also of every local interaction operator hij. For spin 
chains such global states can be generated by using so-called matrix product 
ground states (MPG) J|, |J. A generalization of the MPG concept to arbitrary 
lattices is given by vertex state models [Q— which have been used to construct 
optimum ground states on the hexagonal lattice and on the two-leg ladder. 

In order to construct the global ground state for the present model, we assign 
the following set of vertices to each site on the first sublattice of the Cayley tree: 



era 1 3) 






: IT) 

(5) 

<< ■■ IT) 





X : IT) 



The corresponding vertices on the second sublattice are 




(6) 



Unlike a classical vertex model, each vertex has a single-spin state a\m) as its 
value (or 'weight'), where 

m = — of outgoing arrows — ^ of incoming arrows) . (7) 

The parameters a and a are the same as in (||) . Both sets of vertices differ only 
with respect to the positions of the cr-coefheients. Note that (^|) and are the 
same as on the hexagonal lattice ||, only the global lattice topology is different. 

The global ground state \^>o) is generated by concatenating the vertices at 
all lattice sites. As in usual classical vertex models of statistical physics, the 
connecting bond between adjacent lattice sites is summed out. The generic 
product of vertex weights is replaced by the tensorial product in spin space: 




It can be shown that the resulting global state is indeed an optimum ground 
state of H by collecting all two-spin states which are generated by all possible 
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concatenations of neighbouring vertices^ 
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These 9 two-spin states are orthogonal on all excited local states (H), so (||) 
are the local ground states of /ly. Therefore it is clear that any projection of 
l^o) onto the space of two adjacent lattice sites is a linear combination of local 
ground states. This yields 



for all nearest neighbours i and j and hence also 



0. 



(10) 



(11) 



Since zero is a lower bound of the global ground state energy, the constructed 
vertex state model is indeed an optimum ground state of the global Hamiltonian 
H. 

In contrast to regular lattices there is no canonical way to impose periodic 
boundary conditions on the Cayley tree, so open boundary conditions are used. 
In this case, the vertices on the boundary sites ('leafs') emanate external bonds 
which are not summed out. Independent of the arrow configuration on these 
external bonds, the resulting vertex state model is always an optimum ground 
state of H. Thus the ground state degeneracy grows exponentially with system 
size. 



4 Properties of the ground state 

Each arrow configuration {6} on the external bonds generates an optimum 

ground state l 1 ?^) of H. The calculation of ground state expectation values 
requires taking the average over all these configurations 

< A >*« = 4 E < A >*</> • ( 12 ) 

{b} 

B denotes the total number of such configurations. By using the techniques 
developed in appendices A and B, ground state expectation values of arbitrary 
single-spin observables and two-spin correlation functions can be obtained ex- 
actly. As described in appendix A, the average over all boundary arrow config- 
urations is performed automatically. Note that the formulae given below hold 
for all system sizes, on which the considered observables can be applied (cf. ap- 
pendix B). In particular the results are valid in the thermodynamic limit. This 
pathological effect occurs only on the Cayley tree, optimum ground states on 
regular lattices exhibit a non-trivial finite-size behaviour. 

2 Common prefactors have been omitted. 
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Figure 2: Inverse longitudinal (thin) and transversal (thick) correlation length 
as a function of the parameter a. 



The first interesting expectation values are the components of the canonical 
spin operator 

(Sf)* o = (S!)v o = (S()* o =0. (13) 

This is the expected result as the boundary conditions and the global ground 
states preserve the symmetries of the Hamiltonian, in particular spin-flip sym- 
metry and rotational invariance in the xy-plane of spin space. ( |i~3| ) implies that 
the total magnetization is zero, so the global ground state is antiferromagnetic. 
Although the local magnetization vanishes, its fluctuations are non-trivial 
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(14) 



which increases monotonically from i to | as a function of a 2 , thus covering the 
full range of possible values. Because of rotational symmetry in the a;y-plane of 
spin space we also obtain 

3 3 



((Sffh = ((S?) a >* = ~ 



3 /3 
2 2 +1 



(15) 



In appendix B the transfer matrix technique has been employed to com- 
pute two-point correlation functions. The result for the longitudinal correlation 
function is 



1 + 3a z 
6 + 2a 5 



1 - a' 



(16) 



where \i—j\ is the distance between lattice sites i and j. Note that the correlation 
decays exponentially as a function of the distance and its sign alternates if 
a 2 > 1. The corresponding longitudinal correlation length can be read off from 



equation (n 



In 



Its inverse is given by 
3 + c 



(17) 




Figure 3: Longitudinal nearest-neighbour correlation as a function of the pa- 
rameter a. 



which is plotted as a function of the parameter a in figure |[ Note the divergence 
at a = 1. In this special case, longitudinal correlations between spins with a 
distance of 2 or larger are completely absent (cf. equation (pf)). 

As the longitudinal one, the transversal two-spin correlation function decays 
exponentially as a function of the distance: 

(s-s^f^y.f ^ r'-\ (18) 



Thus the inverse transversal correlation length is 

C*-ta^- (19) 

It is also plotted in figure |2|, together with the longitudinal one. Both inverse 
correlation length, £j — and £ t _1 , are non-zero for all finite values of a, hence the 
model is never critical. 

In the special case = 1 equation ( |l6| ) yields the longitudinal nearest- 

neighbour correlation (SfSf, i)* . Starting at — ^ for a = it decreases mono- 
tonically and approaches — § asymptotically for large values of a, as shown in 
figure [| 

There are three noteworthy special points in the parameter space. The first 
one is a 2 = 1. In this case all non- vanishing vertices of the classical vertex 
model which corresponds to the inner product (^ol^o) (cf. appendix A) have 
the same weight, namely 1. This corresponds to infinitely high temperature in 
the language of classical vertex models, i.e. disorder is maximal. The vanish- 
ing of longitudinal correlations for a 2 — 1 (cf. figure |^) is consistent with this 
interpretation. 

The next interesting special case is the isotropic point a = — V3, a = — 1 
where we can adjust the A-parameters so that the local interaction operator (Q) 
has full SO(3) symmetry. As mentioned in section |[ this model has already 
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been investigated in gj . The reported inverse correlation length £ iso = In 3 for 
open boundary conditions coincides with the results obtained from ( |l7| ) and 

Finally we consider the limit a? — > oo. As can be seen from (§ and (§) the 
global ground state is dominated by only four vertices in this limit, namely 




on the first sublattice and 




on the second one. The two ways, in which these vertices can be combined on 
the Cayley tree, represent the two different Neel states. All spins on the first 
sublattice are in the |3) state, the others are in the |3) state, and vice versa. 
Therefore this special case is called the Neel limit. Note that the ground state 
degeneracy is higher than in the generic case. In addition to the degeneracy due 
to the open boundary conditions there is also a 'bulk degeneracy', as some of 
the local ground states @ become simple tensor products of single-spin states. 



5 Summary 

We have investigated the ground state problem of a class of antiferromagnetic 
spin- 1 models on the Cayley tree with coordination number 3. Apart from the 
lattice topology the Hamiltonian is the same as in our previous works on the 
hexagonal lattice || and the two- leg ladder Q. It is defined in terms of the 
nearest neighbour interaction, which contains 5 continuous parameters and has 
parity invariance, spin-flip invariance, and rotational invariance in the xy-plane 
of spin space. 

Due to the open boundary conditions the ground state degeneracy grows 
exponentially with system size. We have constructed the global ground states 
explicitly by using the vertex state model approach. These are so-called optimum 
ground states, i.e. they are not only ground states of the global Hamiltonian, 
but simultaneously minimize all local interaction operators. The vertex state 
model contains a continuous parameter a, which controls z-axis anisotropy, and 
a discrete parameter a = ±1. 

The calculation of ground state expectation values leads to a classical vertex 
model on the same lattice as the original quantum spin model. Due to the hier- 
archical structure of the Cayley tree this classical model can be solved exactly 
by using recursion relations and the transfer matrix technique. The result of 
our calculations is that the model has vanishing sublattice magnetization and 
exponentially decaying correlation functions. Exact formulae for the nearest 
neighbour correlation, the longitudinal and transversal correlation lengths, and 
for the fluctuations of the magnetization have been derived. For special values 
of a and a the global ground state coincides with the so-called valence bond 
solid (VBS) ground state. 
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Appendix A: Calculation of single-spin expec- 
tation values 



For the moment we consider only one of the many global ground states, namely 
the vertex state model with all boundary arrows pointing out of the leaf sites. 
Denote this ground state by \^q). Within this ground state, the expectation 
value of an observable A il which only acts on the spin at lattice site i, is defined 



as 



{Ai 



(*SI*5> 



(A.l) 



The denominator can be interpreted as two identical vertex state models on 
top of each other, representing the bra- and the ket-vector, respectively. Since 
the vertices at each lattice site generate only local single-spin states, the inner 
product can be taken separately at each lattice site before the interior bonds 
are summed out. Hence (^ol^o) can be interpreted as the partition function of 
a classical vertex model with vertices defined as 



(A.2) 



These vertices have the following properties: 

• The vertex weights are real numbers, not single-spin states. 

• There are two arrow variables on each bond, originating from the bra- and 






1. 



the ket-vector. 

• The vertices are identical on both sublattices of the Cayley tree as a 2 

• Only 20 of the 64 different vertices have a non- vanishing weight since the 
inner product between different S z -eigenstates is zero. 



The numerator of ( A.l ) corresponds to the same classical vertex model as the 
denominator, except for site i where the vertices are modified. At this special 
lattice site, the classical vertices are given by i nser ting the operator Ai between 
the bra- and the ket-vector on the r.h.s. of (A.2). The graphical representa- 
tion of (^ol-^il^o) ls shown in figure ||. The affected site i is surrounded by 
three unmodified branches of the classical vertex model. In the following, these 
branches are denoted by , Z^ . Z^~ , and Z^* , defined as 




(A.3) 



Here the lower index n denotes the order, i.e. the branch contains 2™ — 1 lattice 



sites. Due to the hierarchical structure the values of (A.3) can be calculated by 
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Figure 4: Graphical representation of single-spin expectation values. 



using recursion relations. Each Z n+ \ is given by concatenating two copies of Z n 
to a single classical vertex and summing out the connecting bonds. This yields 
the recursion relations 
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z„ \ + 1Z„ z~ 



zz 



(A.4) 



In order to actually solve these coupled equations it is necessary to specify 
the initial values zf , zf , Z^~ , Z^~* which arc of course determined by the 
boundary conditions. For | \I/q } , which is defined as having all boundary arrows 

pointing out of the leaf sites, the correct choice would be to set zf = 1 and 
the other three initial values to zero^. 

However, in order to calculate ground state expectation values, definition ( |l2| ) 
requires the summation over all arrow configurations on the external bonds. It 
turns out that this summation can be carried out automatically by using the 
initial values 

zf = zf = 1 and (A.5) 
zf = zf=0. (A.6) 



( A.5) ensures that configurations with incoming and outgoing arrows are weighted 



equally and (A.6) eliminates all 'mixed' terms, i.e. terms where bra- and ket- 
vector are different. 

As a consequence of these initial values we obtain 

Zf = Zf and (A.7) 
Zf = zf=0 (A.8) 



3 Multiplying the set of initial values with a common non-zero constant leaves all expecta- 
tion values unchanged. So ZZ^ = 1 is a convenient normalization. 
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for all n > 0. This is immediately clear from the recursion relations. Inserting 

(A.9) 



(|OD, flOD into (|AJ) yields 

Z^n+l — Z^+l = (3 + ft2 ) ( K, 



which has the solution 



zr = v = ( 3 + o ) 



2n(2"-1) 



(A.10) 



Note that the exponent 2 n — 1 is simply the number of lattice sites of the branch. 
It is now straightforward to calculate the ground state expectation value of the 
local observable A as 



(A 



■i/*o 



ZjAj 
Z 



(A.ll) 



Z{Ai) is given by attaching the solution (A. 10) to all bonds of the classical 
vertex modified by the operator (cf. figure^). The denominator is simply 



Z = Z{1) = 2(3 + a 2 ) 



2\N 



(A.12) 



N being the total number of lattice sites. Due to the product structure of 
the solution ( A. 10 ) the results for all expectation values are independent of 



N . Additional powers of 3 + a 2 drop out in the quotient ( A.ll ). As shown in 



appendix B, a similar effect also occurs in the general case of fc-point correlation 
functions. This means that there are no finite size effects. Note that although 
the classical vertex model has this simple product solution, the underlying vertex 
state model generates a highly non-trivial global ground state. 

Equation flA.Sj ) means that the classical vertex model contains no unequal 
arrow pairs on its bonds. The reason for this exact vanishing is that only 
classical vertices with zero or two unequal arrow pairs have a non-vanishing 
weight. Two of the three bonds of a leaf site are external bonds, which do not 
carry unequal arrow pairs (due to the open boundary conditions) . So each leaf 
site provides the next hierarchy of vertices only with equal arrow pairs, and so 
on. Thus unequal arrow pairs can never enter the classical vertex model. This is 
in contrast to the hexagonal lattice || , where unequal arrow pairs are in general 
present, albeit in a very low concentration for most values of the parameter a. 
This is a dynamical effect. 



Appendix B: Calculation of two-spin expecta- 
tion values 



In this appendix we extend the technique developed in appendix A to two-point 



and Bj are observables which only act on 



correlation functions (AiBj)\p . Ai 
lattice sites i and j, respectively. 

The Cayley tree is free of loops, so there is exactly one path from site i to site 
j. Figure g shows the topological structure of the corresponding classical vertex 



model. The branches Z£ v are known from (A.8) and ( A.1C ), so the remaining 
problem is the summation along the path from i to j. If \i — j\ denotes the 
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Figure 5: Graphical representation of two-spin correlation functions. 



distance between these two lattice sites the path is given by the j\ — l)-fold 
product of the periodicity element^] 



(B.13) 



The bonds fiivi and ^2^2 connect this element to the previous and to the next 
periodicity element on the path. The idea is to interpret ( B.13] ) as a transfer 
matrix T [|, ^| and to use its eigenbasis to calculate Tl 1- - 7 ' -1 . 

Each bond carries two binary arrow variables, so T is a 4 x 4-matrix. If the 
following mapping of arrow configurations to matrix indices is used 



=4 



(B.14) 



then the transfer matrix is given by 



1 



T = cz 



V 






2 J 



(B.15) 



cz = (3 + a 2 )^ 2 ™ -1 ) is the prefactor due to the attached branch Z£ v . The 
eigenvalues \k and the corresponding normalized eigenvectors \v,k) of this real 
symmetric matrix are 



X.2 
X3 

Xi 



a )c z 

2\ 



(3 

(1 - a 2 )cz 
2c z 



|ui) = (1, 1,0,0)/V2 

\u 2 ) = (1,-1,0,0)/V2 

|ua) = (0,0,1,0) 

\u 4 ) = (0,0,0,1). 



(B.16) 



This eigensystem can now be used to compute the desired powers of the transfer 
matrix as 



7 1 1»— J I 



E 



x k 



t-il-i 



Uk)(Uk\ 



(B.17) 



4 Strictly speaking the order n of the branch Zf^ can be different for each vertex along the 
path from site i to site j. However, for the same reason as in appendix A the order of these 
branches is completely irrelevant for the calculation of expectation values. 
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The final step is to assemble all components of the classical vertex model as 
shown in figure |. To the left/right-hand side of jB.17p we attach the vertices 



modified by the observable Ai/Bj, respectively. Unmodified branches Z n are 



attached to the remaining four bonds. In analogy to (A.ll) the desired expec- 
tation value is given by 

(AiBj)^ = . (B.18) 

The numerator is the partition function of the modified vertex model calculated 
in this appendix and the denominator is the partition function of the unmodified 
vertex model Z — 2(3 + a 2 ) N . As in the previous appendix, the average over all 
configurations of the boundary arrows is performed automatically. 



Note that the order of the branch Z n in (B.13) enters the calculation only 



via the factor cz which appears in all eigenvalues of T (B.16). However, cz 



drops out in the quotient (B.1S), so the order of the branches attached to the 
path from site i to site j is irrelevant. The same holds for the four branches on 
the l.h.s. and on the r.h.s. of figure ||. Therefore the expectation value (B.18) 



only depends on the observables Ai and Bj and on their distance \i— j\. The 
size of the rest of the Cay ley tree has no influence. 

The technique developed in this appendix can be easily generalized to k- 
point correlation functions. If k is finite there is always a decomposition of the 
full lattice into a finite number of 

• paths, 

• modified vertices, and 

• unmodified branches. 

These objects can be dealt with separately. The resulting expectation values 
depend only on the observables themselves and the length of the paths between 
the vertices modified by these observables. The size of the unmodified branches 
is irrelevant. In this sense there are no finite size effects. 
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